Origin of spin gap in CaV409: effect of frustration and lattice distortion 
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We study the origin of spin-gap in recently discovered ma- 
terial CaV409. We analyze the spin-1/2 Heisenberg model 
on the 1/5 depleted square lattice with nearest neighbor (nn) 
and next nearest neighbor (nun) interactions, in terms of the 
singlet and triplet states of the 4-spin plaquettes and 2-spin 
dimers. Phase diagram of the model is obtained within a 
linear "spin- wave" -like approximation, and is shown to agree 
well with the earlier results of QMC simulations for nn inter- 
actions. We further propose that the special lattice structure 
of CaV409 naturally leads to lattice distortions, which en- 
hances the spin-gap via a spin-Peierls mechanism. 

PACS: 75.10.Jm, 75.30.Kz, 75.40.Cx, 75.50.Ee 



Recent discovery Q of a quantum disordered phase and 
spin-gap in the layered magnet CaV409 has attracted 
considerable interest ||,^;^,||j6j. The magnetic system 
can be described by a Heisenberg model for spins of Vana- 
dium ions {S = 1/2) on a 1/5-depleted square lattice. At 
each site of this bipartite lattice three bonds meet: two of 
them belong to the 4-spin plaquettes covering the lattice 
(plaquette bonds), whereas the third one (dimer bond) 
connects a plaquette with its neighbor (see Fig. 1). Since 
the coupling between spins is mediated by superexchange 
via intermediate Oxygens, a strong next nearest neighbor 
interaction is also expected ^j^. 

We are thus led to the following Hamiltonian: 
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where the nn interaction J„„ equals Jg ( Ji) for plaquette 
(dimer) bond. It is evident that this model has disordered 
singlet ground states in two limits: for Jg 3> Ji,J2 the 
ground state is a product of singlets on each plaquette, 
and for Ji ^ Jq, J2 it consists of singlets on the dimers. 
However, a physically relevant choice of exchange param- 
eters is Jq M Ji. 

For model (^ a spin-gap can arise either from the non- 
equivalence of the dimer and plaquette bonds (Jo ^ Ji), 
or from the effects of frustration due to nnn couplings. 
Previous theoretical works have mainly focused on the 
former. Quantum Monte Carlo (QMC) simulations, done 
with J2 = show that the 1/5-depleted square lattice 
has Neel order for Jq = Ji. However, the system is close 
to the transition into the disordered phase. 

The aim of this paper is twofold. First, we describe an 
efficient analytical approach to study the T = phase 



diagram of (|l|) in the Ji/ Jq — J2/J0 plane, and show that 
a moderate J2 ^ 0.2Jo, with Ji = Jq, can account for 
the experimental data. Second, we point out that the 
special lattice structure of this system inevitably leads 
to lattice distortions, resulting in the strengthening of 
the Plaquette bonds with respect to the Dimers. Thus 
the spin-Peierls mechanism cooperates with the intrinsic 
tendencies of the system in forming a spin-gap. 

The phase diagram of the Hamiltonian (|l|) for J2 = 
has been studied analytically using linear spin- wave 
strong coupling perturbation [^|J^ , mean-field Schwinger 
boson and cluster |H approaches. As expected, both 
spin-wave and Schwinger boson theories substantially 
overestimate the region of stability of the Neel phase com- 
pared to QMC data |^ , while expansions around strong- 
coupling limits favor disordered phases. Here, we inves- 
tigate this model using bosonization techniques which 
correctly account for short-range spin correlations inside 
the plaquette and dimer blocks. 

For the closely related case of the frustrated square- 
lattice antiferromagnet the disordered spin-state is a 
dimer state, formation of which is accompanied by a 
spontaneous breaking of the lattice symmetry In 
contrast, for the CaV409-lattice the choice of spins which 
form singlets is determined by the structure of the lat- 
tice and there is no spontaneous symmetry breaking. 
However, two different types of disordered short-range 
RVB states, with spin-singlets formed on plaquettes and 
dimers, are possible |^,^. For different values of model 
parameters one needs to consider representations for spin 
operators in terms of both dimer and plaquette states. 
We generalize previous derivations of such representa- 
tions @,|,|lO|,0 for the two cases. 

The starting point of these representations are non- 
interacting spin blocks. Let states of a single block 
be given by \a). In case of dimers, they are a singlet 
\s) {Eg — ~jJi) and a triplet \ta), a = x,y, and z 
[Et = \ Ji)- AH 16 states of a four spin plaquette can be 
found in Ref. [||J|]. The lowest levels, once again, are a 
singlet with energy _Es = — 2Jo -I- 5 J2, and a triplet with 
Et — — Jq ^ \ J2- In the bosonic representation for pla- 
quette spins we will restrict ourselves to these four states 
only. This assumes that occupation numbers of all higher 
levels are small. 

The site spins are expressed in terms of the basis 
block states as 
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s, = («|s,|/?)z"^ (2) 

where Z°'^ is the projection operator |a)(/3| and summa- 
tion over repeated indices is assumed. 

Consider first the matrix elements that occur in Eq. , 
in the subspace of one singlet and three triplet states. Us- 
ing rotational invariance in spin space and time-reversal 
symmetry one gets 

(s|Sr|s)=0, {s\St\tp) = 5^pA\^, 
(i„|Sf|t^)=*e"^^^j„ (3) 

where e"^^^ is the totally antisymmetric tensor and A\^^ 
are real. Let each spin of the block be decomposed 
as Sj = Lj + Mi, where Lj has nonzero matrix ele- 
ments between triplets and singlet, while acts be- 
tween triplets only. It is easy to see that Mj must be 
proportional to the total spin of the block, and, thus, 
is independent of i, whereas cx (—1)* for blocks con- 
sisting of equivalent spins. Further calculations using 
explicit forms of singlet and triplet states are straightfor- 
ward and give A\^ = (— 1)Y2, A\i = 1/2 for dimers, and 
Al^ = (-l)V%/6, A\^ ^ 1/4 for plaquettes. 

We define the vacuum |0) and four boson operators 
that yield the four physical states \s) = s^|0), \ta) — 
t+|0). The projection operators are naturally expressed 
as Z'^^" — s^ta, Z*"*" = t'^tp, and so on. Block spins 
represented via these boson operators are 

5f = ^{s+t^ + t+s) ~ ^e'^^^t+t^ , for dimers, 

5*^ = i-^(s+<„ + t+s) - \e"'^''t+t^ , for plaquettes. 

Commutation relations between spins are satisfied as 
long as bosonic representation preserves the algebra of 
the projection operators, i.e. the relation Z"^ Z" ^ — 
Sa'pZ"^ . This requirement restricts the number of 
bosons allowed on each block to one: 

s+s + t+tc, - 1 . (5) 

With the help of this constraint the Hamiltonian of a sin- 
gle block becomes Hb = EgS+s + Ett+ta- There are two 
ways to implement the constraint (H). First is the slave- 
boson treatment of Sachdev and Bhatt , where relation 
is imposed by a Lagrange multiplyer. The second ap- 
proach [0 is to express the s-operators in terms of the 
ta-operators 

s+ = s = y/l - t+to, , (6) 

and substitute into Eq. (|^). Instead of the Holstein- 
Primakoff type representation (^), we can also consider 
an analog of the nonhermitian Dyson-Maleev representa- 
tion: s+ = 1 and s = 1 - t+ta Note that the block 
Hamiltonian is treated identically in both approaches. 
Differences appear when interactions between the blocks 
are switched on. 



As in case of spin-wave theory for ordered magnetic 
phases, one might expect even the linear approximation, 
which neglects interaction between excitations, to work 
well. This approximation can be formulated for both 
types of block spin representations in the same way. It 
consists of replacing s and s+ by 1, when calculating 
(Si • Sj) for pairs of spins from different blocks. Also, 
only terms of second order in triplet operators should be 
kept. Diagonalization of the resulting quadratic form is 
done by a standard Bogoliubov transformation, and one 
finds a 3-fold degenerate spectrum of triplet excitations. 
For small coupling between the blocks the spectrum is 
positive with a gap. Increasing the interaction between 
the blocks decreases the gap, which finally vanishes at 
the transition between disordered and ordered phases. 

First, consider the plaquette singlet phase which ex- 
ists for large Jq. The spectrum of spin-1 excitations is 
threefold degenerate and has the dispersion 

w^(k) = Jo [Jo + - 2J2)(cos K + cos ky)] . (7) 

The minimum of the spectrum is at (tt, tt) for (Ji — 2J2) > 
and at (0, 0) for ( Ji - 2 J2) < 0. From Eq. (0) one finds 
the region of stability of the plaquette phase, shown in 
Fig. 2. At J2 = 0, singlets on 4-spin plaquettes become 
unstable at the critical ratio Jo/Ji|cr — |, which is not 
too far from the QMC estimate Jo/Ji|cr — 1-1 @- The 
total energy of this phase, per spin, is 

Els. - -kJo + lJ2 + i^Y. ~ M ■ (8) 

k 

It consists of the energy of noninteracting plaquette sin- 
glets and the energy of zero-point fluctuations. 

In the dimer state each crystal unit cell has two dimers. 
Therefore, there are two different branches of S" = 1 
magnons in the Brillouin zone. However, calculations are 
greatly simplified if instead we consider only one type of 
dimers, which are defined in the new Brilloin zone corre- 
sponding to the lattice formed by the centers of dimer 
bonds. This procedure is quite similar to performing 
a spin-wave expansion around the rotating quantization 
axis, as is often done in standard spin- wave theory. As a 
result, we obtain one triply degenerate excitation mode 
in the new Brillouin zone, which is twice the original one, 

ojjCk} = Ji [Ji - (Jo - J2){costj, - cosfcj^) 

- J2 cos{ka: + ky)] . (9) 

The minimum of the spectrum is at k = (0, tt) for 
J2 < ^^0, and moves into the interior of the zone for 
larger J2. At J2 = 0, the dimer phase is unstable for 
Jo/Ji > while the corresponding critical ratio from 
QMC is Jo/Ji|cr = 0.6 ± 0.05 @. Note that the insta- 
bilities of the Neel phase at J2 = 0, as predicted by our 
linear bosonic theory, coincide with the mean field clus- 
ter approach B. However, the latter is not appropriate 
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to study spin-liquid phases. The total energy per spin in 
the dimer phase is 



El 
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(10) 



Phase diagram of the model (|l|) follows from Eqs. (|^) 
and (||). When frustration exceeds the critical value 
J" = O.25J0, the Neel phase ceases to exist. A first or- 
der transition line separates plaquette and dimer phases 
since their symmetries are different Q . This line is found 
by comparing (||) and ( p^ and is shown in Fig. 2 [ p^ . 

As J2 grows, the energies of the omitted plaquette 
states decrease making the linear approximation less sat- 
isfactory, and at J2 = Jo the second singlet which consists 
of two crossing dimers becomes the ground state of the 
4-spin plaquette We have checked that this phase 

is not stable in the linear approximation for any values 
of the parameters. Another possible short-range RVB 
state is the plaquette-RVB on large squares formed by J2 
bonds, but it was also found to be unstable. Therefore, 
for large values of J2 magnetic order should be stabilized 
again. It is easy to see that in the limit J2 S> Joi Ji spins 
are arranged into two interpenetrating Neel ordered sub- 
lattices which are decoupled at the classical level. The de- 
generacy with respect to a relative orientation of antifer- 
romagnetic vectors should be removed by quantum fluc- 
tuations, providing an example of 'order-from-disorder' 
phenomenon |13|. We have presented in Fig. 2 transition 
line between this ordered state and disordered plaquette 
phase. 

A theoretical comment is in order here. Calculations 
described above can be repeated within an alternative 
slave-boson treatment 0, which amounts to introduc- 
ing a singlet condensate 7^ 1 and a chemical poten- 
tial ^, and finding their values self-consistently. Solv- 
ing the self-consistent equations for the phase bound- 
aries, we found, at J2 = 0, that dimer phase is stable 
for Jq/Ji < 0.874, whereas plaquette-RVB exists up to 
Jo/Ji = 0.76. Thus, unlike linear approximation, mean- 
field slave-boson treatment predicts no ordered phase 
even at J2 = 0, in striking contradiction with numer- 
ical simulations and linear approximation results. Sit- 
uation is not improved even upon taking triplet-triplet 
interaction into account - its effect on the location of 
phase boundaries is ridiculously small, less then 1 per 
cent. These findings clearly favor our linear approxima- 
tion as the most suitable one for the problem at hand. 

To make contact with experiments, let us first confine 
ourselves to Jq = Ji. Note that only second terms of 
Eq. (Q) contribute to the uniform susceptibility, leading 
to X ~ ^j^n(a;k)(l + "(^k)), where n{ui) is the Bose 
factor. We can compare the Curie- Weiss relation 9 — 
I (Jo + J2) and the gap A, in Eq. appropriate for the 
PRVB phase, with the experimental Curie Weiss constant 
and the gap determined from the T-dependence of the 



susceptibility. We find that the experimental numbers 
A = 107 K and 6 = 220 K ||l| can arise from two sets of 
exchange constants (i) Jo = 245 K, J2 — 48 K and (ii) 
Jo = 193 K, J2 = 100 K. Thus this theory, by itself, 
can provide a consistent description of the properties of 
CaV40g observed so far. 

However, we would now like to point out an additional 
aspect of this material which has so far gone unnoticed. 
Because of the special lattice structure, spin-phonon cou- 
plings will cause a gain in magnetic energy which is lin- 
ear in the lattice distortion, whereas the loss in elastic 
energy is always quadratic. Thus there will always be a 
lattice distortion in this system which will cause the pla- 
quette bonds to shrink and the dimer bonds to elongate. 
This will lead to Jq > Ji. This spin-Peierls mechanism 
will enhance the stability of the PRVB phase and con- 
tribute to the spin-gap. 

To get a quantitative estimate, let the shrinking of 
the plaquettes change the distance between the Vana- 
dium ions from R to R ± 2u. Using the phenomeno- 
logical relation J{R) ~ const /R^'^ the exchange 

constants Jq (Ji) increase (decrease) by 5J = 20uJ/R. 
The magnetic energy gain per Vanadium atom is Cm = 
—SJ{Ci — O.5C2), where Ci and C2 are spin correlations 
on the plaquette and dimer bonds in the absence of dis- 
tortion. It is sufficient to evaluate them on the Neel state, 
where they both are —1/4. Thus the gain in magnetic 
energy per unit volume is Em ~ {20uJ/8R){l/R'^R±), 
where R± is the distance between magnetic layers. A 
rough estimate for the elastic energy per unit volume is 
Sph ~ \B(u/Rf, where B is a bulk modulus of the 
material. Minimizing the energies leads to a distortion 
u ^ {5J/2BRR±), which induces a variation of the ex- 
change integrals by 



SJ/Jfn [hOJ/R^R^B). 



(11) 



Taking i? 3 A, i?^ = 5 A, J = 300 K, B = 
10^2 dynes/cm^, we get u/R - lO'^ and 5J/J ~ lO'^- 
10^^. These values are similar to what is found in most of 
the spin-Peierls materials (l^, where Jq/Ji ^ 1.3. This 
similarity is not unexpected as the linear gain in mag- 
netic energy is somewhat analogous to u^f^ gain in a ID 
spin-i chain Jl6t . 

If J2 was negligible, and the spin-Peirels mechanism 
was entirely responsible for the spin-gap, we would use 
Eq. (^ with J2 = and the Curie- Weiss relation 6 = 
■|(2Jo -l- Ji) to compare with experiments. This leads to 
estimates Jo = 331 K, Ji = 218 K. The ratio Jq/Ji ^ 
corresponds to 5 J — 0.2 J, which is somewhat larger than 
our earlier back of the envelop estimate, suggesting that 
some J2 is needed to explain the data. 

We note that here the spin-Peierls mechanism does 
not lead to a spontaneous breaking of the lattice sym- 
metry. Hence, we do not expect a sharp change in 
the lattice distortion u with temperature. Instead, it 
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should follow short-range order and the spin-gap should 
be temperature-dependent. Finally, we point out that 
most known spin-Peierls materials have first order struc- 
tural phase transition at high temperature jist - This 
creates a soft phonon mode in the system, which in the 
end favors the spin-Peierls phenomena [0. Curiously 
enough, a small discontinuous jump in xC^) ^-t 340 K 
was observed also in CaV409 but was attributed to 
an admixture of the VO2 phase. In view of our proposal 
this point should be studied more carefully. 

Our discussion of possible spin-Peierls phenomena is 
only an order of magnitude estimate. Nevertheless it 
shows that spin-phonon coupling may successfully "co- 
operate" with intrinsic tendencies for forming a spin-gap 
due to the frustrating next-nearest-neighbor interaction 
present in CaV409. 
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FIG. 1. Lattice structure of CaV409. Three types of ex- 
change bonds are indicated by thick lines. The pattern of 
lattice distortion is shown schematically by thin dashed lines. 

FIG. 2. Phase diagram of the model in linear approxima- 
tion. Thick (thin) solid line denotes second (first) order phase 
transition. Regions of stability of dimer (plaquette) phase are 
shown by long (short) dashed lines. 
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FIG. 1 of Starykhetal. 




FIG.2 of Starykh et al. 



